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ABSTRACT 

O 

o 

CN ' Five years of Rossi X-ray Timing Explorer (RXTE) observations of the 

Oh' Galactic black-hole candidates IE 1740.7-2942 and GRS 1758-258 show a pe- 

^ ■ riodic modulation with amplitude 3-4% in each source at 12.73 ± 0.05 dy and 

• 18.45 ±0.10 dy, respectively. We interpret the modulations as orbital, suggesting 

that the objects have red-giant companions. Combining the RXTE data with 
earlier data (Zhang et al. 1997) from the Burst and Transient Source Experiment 
on the Compton Gamma-Ray Observatory, we find a long period or quasi-period 



> 

in 



'n . of about 600 dy in IE 1740.7-2942, and a suggestion of a similar 600-dy period 

Q"^ ■ in GRS 1758-258. These timescales are longer than any yet found for either pre- 

^ ■ cessing systems like Her X-1 and SS 433 or binaries like LMC X-3 and Cyg X-1 

O ■ with more irregular long periods. 

Subject headings: x-rays:stars — stars, individual: (IE 1740.7-2942) — stars,individual:(GRS 
S ■ 1758-258) 



1. Introduction 

The Galactic-bulge x-ray sources IE 1740.7-2942 and GRS 1758-258 are generally 
called black-hole candidates due to the similarity of their x-ray spectral and timing be- 
havior to that of Cyg X-1 in its usual hard state. Like Cyg X-1, both sources occasionally 
enter an intermediate or soft state, but the evolution of their spectral hardness and luminos- 
ity is very different (Smith, Heindl & Swank 2002). Both have prominent, bright radio lobes 
about an arcminute in size (Mirabel et al. 1992; Rodriguez et al. 1992), while Cyg X-1 does 
not (Marti et al. 1996), showing only a milliarcsecond jet near its core (Stirling et al. 2001). 
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The counterparts of IE 1740.7-2942 and GRS 1758-258 are unknown due to high ex- 
tinction, therefore there are no orbital solutions or estimated masses. For GRS 1758-258, 
Heindl & Smith (2002) recently used Chandra data to confirm the association of the x- 
ray source with the radio core source ("VLA-C") and Cui et al. (2001) did the same 
for IE 1740.7-2942. Marti et al. (1998) identified two candidate counterparts to VLA- 
C/ GRS 1758-258 in /- and K-hand images. The brighter and closer candidate was found 
through multi-band photometry and near-infrared spectroscopy to be a likely KO III giant. 
Revised astrometry (Rothstein et al. 2002) of infrared observations by Eikenberry et al. 
(2001a) confirm that this star ("star A") is consistent with VLA-C at the 3(t level. Marti et 
al. (2000) and Eikenberry et al. (2001a) agree on several possible high-mass candidates for 
the companion of IE 1740.7-2942 in its more crowded and obsciued field, but a low-mass 
companion would be unobservable at the current K-hand sensitivity. 



2. Observations 

We use a five-year (1997-2001) series of observations by the Rossi X-ray Timing Explorer 
(RXTE) of IE 1740.7-2942 and GRS 1758-258 (Main et al. 1999; Smith et al. 2001; Smith, 
Heindl & Swank 2002). The observations, each of 1000-1500 s, were taken approximately 
weekly in 1997-2000 and twice weekly since early 2001. Monthly observations in 1996 are 
too sparse to improve our results and are not included. Data cannot be taken from late 
November to late January of each year when the Sun is close to the Galactic Center. 

We use the Proportional Counter Array (PCA), layer 1, in the range 2.5 to 25 keV. 
Instrumental background has been subtracted using the "faint source" model and Galactic 
diffuse emission has been subtracted using pointings to nearby fields without bright point 
sources. To compensate for gain changes, we accumulate counts in bands of constant energy, 
not channel number. We offset-point by about a half degree from each source to avoid 
nearby bright sources. Details of the offsets and background pointings are given in Main et 
al. (1999). 

Since our orbital signals are small, we made two checks for systematic errors. Different 
observations use different subsets of the 5 detectors of the PCA. We repeated the entire 
analysis below using only the third detector, which is always on, and found no change except 
for the expected statistical degradation. Shght changes in the PCA field of view due to 
changes in the roll angle of the spacecraft take place over months, and so cannot mimic the 
shorter orbital periods. 
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Fig. 1. — Analysis of the lightcurve of IE 1740.7-2942. a): Count rate > 12 kcV per detector 
of the PCA. The error bars are smaller than the plotting symbols. The smoothed curve is 
described in § 3. b): The same data high-pass filtered by subtraction of the smoothed curve, 
c): Fine histogram: Lomb-Scargle periodogram of the data in panel b). Coarse histogram: 
highest values achieved by applying the same analysis to scrambled data (see text) . d) : Data 
from panel b) folded on the period of the peak in panel c), and normalized to the average 
count rate. The smooth curve is the best fit sine function. 
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Fig. 2. — Same as Figure 1 for GRS 1758-258. GRS 1758-258 spent more of its time in the 
soft state, so more data have been removed. 
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3. Analysis and results 

3.1. Orbital modulation 

During long intervals of relatively stable hard state emission, these sources show low- 
amplitude modulations that are most naturally interpreted as orbital. Figures la and 2a 
show the 1997-2001 background-subtracted lightcurves from 12-25 keV. There are gaps 
due both to the annual solar constraint and to transitions to the intermediate and soft 
states (Smith et al. 2001; Smith, Hcindl & Swank 2002). Luminosity variations are so large 
during the state transitions and soft periods that they decrease the significance of the orbital 
measurements. The data cut is impartial in that only the spectral power-law photon index 
is used to determine which parts of the data to remove: it is harder than 2.0 in the surviving 
data. 

We high-pass-filtered the data by subtracting out a smoothed version of the lightcurve 
(Figures la and 2a) leaving only high-frequency residuals (Figures lb and 2b). To make the 
smoothed curve, we replaced each data point with a value generated by fitting a polynomial 
of order P to all data within TV days of it. Raising N or lowering P increases the amount 
of smoothing. The figure shown is for N — 1{) and P — 1, but the results are not highly 
sensitive to these values. 

We then took a Lomb-Scargle periodogram of the residuals with the result shown in 
Figures Ic and 2c: peaks at 12.73 ± 0.05 dy in IE 1740.7-2942 and 18.45 ± 0.10 dy in 
GRS 1758-258. The errors are the half width at half maximum of the peaks. Although 
these values would be near the Nyquist frequency for weekly sampling, the irregularity in 
the sampling times extends the useful frequency range many times higher. 

Figures Id and 2d show the residuals folded on the best-fit periods and divided by the 
average count rate. The phases are referenced to Oh UT on 1 January 1996. The best-fit 
amphtudes are (3.43 ± 0.26)% for IE 1740.7-2942 and (4.42 ± 0.32)% for GRS 1758-258. 

When the data from each source are divided into three roughly equal intervals, the 
orbital modulation is seen in all three as the highest periodogram peak in the range 1-30 dy. 
We also repeated the analysis in 2.5-4.0 keV, 4.0-6.0 keV, 6.0-8.0 keV, 8.0-10.0 keV, and 
10.0-12.0 kcV bands. The modulations are not as strong, but in every band they are the 
highest peaks from 1-30 dy. Neither the amplitude nor phase varies significantly with energy, 
although amplitude may increase slightly with energy in IE 1740.7-2942 (Figure 3). The 
2.5-4.0 keV band is not used for IE 1740.7-2942 due to its high absorption column. 

We demonstrated by simulations that our process neither creates spurious signals nor 
destroys real ones. To look for spurious signals, we randomly permuted the values of each 
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Fig. 3. — Dependence of the parameters of the best-fit orbital sinusoids on energy. Top: 
phase (at Oh UT on 1 January 1996). Bottom: amphtude. 
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triplet of data points in Figures la and 2a. This creates data with the same samphng 
times, average, and long-term evolution as the real set but with any high-frequency signals 
destroyed. A thousand differently-permuted datasets were analyzed for each source, giving 
the coarse histograms in Figures Ic and 2c, which show the highest power seen in one-day- 
wide period bins in any of the 1000 trials. It is clear that peaks the size of the real signals 
do not occur. 

To test the usefulness of the high-pass filtering, we added to the original data small 
sinusoidal signals with 3% to 20% of the mean flux and 10-20 dy periods. Wc then did 
the analysis with and without filtering. We found that filtering improves the sensitivity 
to artificial signals by a factor of 2 to 4, varying with the period and data set. Artificial 
signals with the amplitude of the real ones can be seen, but only by filtering first, proving 
the necessity of the procedure. 

3.2. Super-orbital period modulation 

To search for longer periods, we combined data (Zhang et al. 1997) from the Burst and 
Transient Source Experiment (BATSE) on the Compton Gamma-Ray Observatory with our 
RXTE/PCA data, extrapolating the PCA data to the BATSE range of 20-100 keV. For 
most black holes in the hard state this is not a bad approximation: only beyond 100 keV do 
the spectra steepen significantly. Figures 4 and 5 show the combined lightcurves and their 
periodograms. The major dips in the lightcurves all correspond to softenings of the spectra. 

The periodogram for IE 1740.7-2942 shows a clear peak around 600 dy. Each individual 
600-dy cycle is visible if the single low point from BATSE in 1993 is included. The data for 
GRS 1758-258 are less conclusive. The BATSE data are poorer for GRS 1758-258, due in 
part to its slightly lower luminosity but probably also to the proximity of the bright Z-source 
GX 5-1. 



4. Interpretation 

If the 12- and 18-dy modulations are orbital, then the companions of IE 1740.7-2942 
and GRS 1758-258 must be giant stars to fill their Roche lobes. For GRS 1758-258, infrared 
observations rule out a high-mass companion (Marti et al. 1998; Eikenberry et al. 2001a) 
and therefore make wind accretion less likely. Star A of Marti et al. (1998), a KO III giant, 
would almost exactly fill its Roche lobe with a IOMq companion and an 18-dy orbit (Smith, 
Heindl & Swank 2000; Rothstein et al. 2002). Perhaps both systems have companions on 
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Fig. 4.— BATSE fluxes for IE 1740.7-2942 from Zhang et al. (1997) (crosses) and PCA 
data extrapolated to the same energy range (diamonds). Top: flux vs. time, 20-100 keV. 
Bottom: Lomb-Scargle periodogram. 
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Fig. 5.— Joint BATSE/PCA lightcurve and periodogram for GRS 1758-258. 
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the first giant branch and are in a long-lived evolutionary phase with accretion driven by the 
gradual nuclear evolution and expansion of the secondaries, as has been suggested for long- 
period neutron-star binaries like Cyg X-2 (Webbink, Rappaport & Savonije 1983). King 
et al. (1997), however, found that such systems should have unstable disks and appear as 
transients if the primary is a black hole. 

Low-amplitude (~7%) x-ray modulation with a similar shape was recently discovered 
in another persistently luminous black-hole candidate, LMC X-3, at its previously-known 
orbital period of 1.7 dy (Boyd, Smale & Dolan 2001). Cyg X-1 also shows low-amplitude 
orbital x-ray modulation, which is not observed in the soft state. Wen ct al. (1999) offered 
two interpretations: partial absorption by the optically thin wind of the secondary (sec also 
Brocksopp et al. 1999), and partial obscuration of an extended emission region (corona) by 
the accretion stream. The lack of modulation in the soft state was explained in the first 
picture by increased ionization and decreased mass in the wind, and in the second picture 
by shrinkage of the corona. 

Because their amplitude and phase are independent of energy (Figure 3), the modula- 
tions of IE 1740.7-2942 and GRS 1758-258 are consistent with an extended emission region 
and an optically thick absorber such as the companion star. Wind absorption would produce 

more modulation at lower energy, as is seen in Cyg X-1 (Brocksopp et al. 1999). Assuming 
the secondary itself is the occulter, then given the size and inclination of the binaries we 
could determine the size of the hard-state emission region from the amplitude of the mod- 
ulation. We arc planning infrared observations of star A over the orbit to determine the 
primary mass, system size, and inclination using radial velocity, velocity broadening, and 
ellipsoidal modulation measurements. 

The prominent radio lobes of IE 1740.7-2942 and GRS 1758-258, appearing far from 
the core sources, imply strong collimation of the jets. Spruit, Foglizzo, & Stehle (1997) show 
that in a jet coUimated by the poloidal magnetic field of an accretion disk, collimation is 
proportional to the ratio of the outer to inner disk radii. The very large disks implied by 
the long binary periods in these systems may support this model. 

The 600-dy super-orbital periods we report here arc the longest seen in any system 
so far. Her X-1 and LMC X-4 have very regular long x-ray periods, about 35 dy and 
30 dy respectively (Tananbaum et al. 1972; Lang et al. 1981), which are explained by disk 
precession with periodic partial obscuration of the x-ray emitting regions (e.g. Wijers & 
Pringle 1999). Jet precession at a period of about 162 dy is observed via Doppler shifts 
of emission lines in SS 433 (Margon et al. 1979; Eikenberry et al. 2001b), and is presumed 
to be caused by precession of the disk. Other systems, such as the black-hole candidate 
LMC X-3, show more irregular variations, associated with spectral state changes that are 



- 11 - 



not consistent with disk obscuration (Paul, Kitamoto & Makino 2000; Wilms et al. 2001). 
Brocksopp, Groot, & Wilms (2001) review recent models of long-period variability involving 
changes in accretion rate either due to a limit cycle caused by x-ray evaporation of the outer 
parts of the accretion disk (Shields et al. 1986) or changes in the companion which cause it 
to fill its Roche lobe only occasionally (Wu et al. 2001). Cyg X-1 in different epochs has 
been reported to show long periods at 294 dy and 142 dy (Priedhorsky, Terrell, & Holt 1983; 
Brocksopp et al. 1999). 

Ogilvie & Dubus (2001) calculated the susceptibility of disks in x-ray binaries to radiation- 
driven warping. We find the binary separations of IE 1740.7-2942 and GRS 1758-258 to 
be about (2.5-5) x 10^^ cm for black hole masses of 3-20 Mq. Using Figure 7 of Ogilvie & 
Dubus (2001), we find that for a secondary of 1.1 Mq (i.e. star A), black hole masses greater 
than 6 Mq fall in the "indeterminate instability zone" , with possible cycling between flat and 
warped disks. Such cycling might explain long, semi-regular periods with spectral changes, 
even with constant mass transfer from the companion. For lower masses of the primary, the 
systems would fall in the regime of persistent warping and stable precession, as in SS 433. 

A recent observation of GRS 1758-258, however, suggests that variations in accretion do 
play a role in its long-term variability. An extreme luminosity drop is visible at about 2002.25 
in Figure 5. Smith et al. (2001) found that GRS 1758-258 entered a soft state suddenly and 
then decayed exponentially with a one-month timescale. Exponentially decaying emission in 
the soft state is easily explained by the steady draining of a thin accretion disk with little 
or no mass input. 

This work was supported by NASA grant NAG5-7265. 
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